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An exact solution is presented for the frequency-dependent charge susceptibility of the spinless 
Falicov-Kimball model by using dynamical mean-field theory. We develop a nontrivial application 
of the Baym-Kadanoff "conserving approximation" formalism to exactly determine the frequency- 
dependent vertex function (which turns out to assume a particularly simple form) . We show how the 
static and dynamic susceptibilities are decoupled in this model and how the dynamic susceptibility 
generically does not show any signal of the low-temperature charge-density-wave phase transition. 
We also examine the temperature evolution of the dynamic charge susceptibility for the special case 
of half-filling. 



I. INTRODUCTION 



The Falicov-Kimball modefl was introduced in 1969 as a thermodynamic model for metal-insulator transitions in 
systems that have two different types of electrons — itinerant conduction electrons and localized /-electrons. Twenty 
years later, Brandt and Mielscho solved the Falicov-Kimball model exactly in the limit of infinite dimensions by 
using dynamical mean-field theory and the Baym-Kadanoff formalism. They also showed how to calculate static 
susceptibilities and found an fsing-like phase transition to a chessboard charge-density-wave phase at half-filling. 
Freericksu then showed that the system also displayed incommensurate order and phase separation at other fillings. 

The Hamiltonian of the spinless Falicov-Kimball modetl is 
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where d\ (di) is the spinless conduction electron creation (annihilation) operator at lattice site i and Wi = or 1 is 
a classical variable corresponding to the localized /-electron number at site i. The hopping matrix between nearest 
neighbors (on a hypercubic lattice in d-dimensions, with d — » oo) is — t* with t* chosen as our energy 

unit, e/ is the localized electron level, fj, is the chemical potential and U is the mutual electron repulsion when a 
conduction electron and a localized /-electron both occupy the same lattice site. 

In this contribution we examine the d-electron dynamical charge susceptibility [where nj(r) = 
exp(ri?)ni(0) exp(— tH) and rii = d\di\ 
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and its analytic continuation to the real-frequency axis. Here ivi = 2mTl is the Bosonic Matsubara frequency, Z is 
the partition function, and Rj is the position operator for lattice site j. The Falicov-Kimball model is the simplest 
many-body problem that has nontrivial dynamics. The model is simple enough that an exact solution can be found 
for the charge susceptibility, but is complicated enough to show many-body effects. Shvaikaa has shown how to 
determine the lattice charge susceptibility from a diagrammatic analysis of the atomic problem. Here we provide an 
alternate derivation, stressing the Baym-Kadanoff approach (which produces a much simpler form for the vertex), 
and we provide numerical results to complement the formalism. 

We begin with a discussion of some special symmetries of the Falicov-Kimball model. Since the total conduction 
electron number is a conserved quantity, we find that Y) j Tij(r) — Y]j Tij(0) has no r-dependence. Hence, the uniform 
charge susceptibility (q=0) vanishes for all nonzero frequencies. Similarly, the local /-occupation is also conserved 
[wi(r) = uii(0)], since Wi commutes with the lattice Hamiltonian. This implies that both the //-charge susceptibility 
and the mixed df-charge susceptibility have no r-dependence and therefore are static, with no frequency dependence. 
Furthermore, a simple application of the chain rule and the definition of each of the susceptibilities as a derivative of 
the corresponding average electron filling with respect to an external field, shows that there is a discontinuity in the 
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dd-charge susceptibility at zero frequency^ due to the coupling of the d and /-electrons. This implies that there is a 
decoupling between the static and dynamic dd-charge susceptibility! |-,r-,|-| 

The solution of the Falicov-Kimball model has been outlined in detail elsewhere&Ert]. Here we summarize the main 
points to establish our notation. The local Green's function at the Fermionic Matsubara frequency iuj n = inT(2n + 1) 
is defined by 



P . I e -PH, 



'd(r)d+ (0)5(A)) 



G n = G{iu n ) = -Tr d Tr t T T J dre w - T ^ v ^ w w/ , (3) 

with 

Z=Z (n)+e-^f-^Z (^-U), (4) 

the atomic partition function expressed in terms of 

ZoOt) = Tr d (e- /3 ^5(A)), H Q = -pePd, (5) 

In the above equations, the atomic Hamiltonian H a t is the Hamiltonian of Eq. ([!]) restricted to one site, with t* = 0, 
and all time dependence is with respect to this atomic Hamiltonian. The evolution operator S(X) satisfies 
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Jo Jo 



(6) 



with A(r — t') a time-dependent atomic field adjusted to make the atomic Green's function equal to the local lattice 
Green's function. We define an effective medium by 

GQ l (iw n ) = G^ 1 + £„ = iuj n +fi-X„, (7) 

with £„ the local self-energy and \ n the Fourier transform of A(r). The trace in Eq. (^) can be evaluated directly to 
yield 

Gn = woGo{iuj n ) + w±[Ga l {iu n ) - U}' 1 , (8) 

with wq = 1 — w\ and 

wi=exp[-P(e f -n)]Z (v-U)/Z. (9) 

The self-consistency relation needed to determine A„ and hence G n is to equate the local lattice Green's function to 
the atomic Green's function via 



G n = I de P(£) , (10) 

J_ 00 iu) n + [i - E n - e 

with p(e) = exp(— 6 2 )/y / 7t the noninteracting density of states for the infinite-dimensional hypercubic lattice. 

The iterative algorithm to solve for G n starts with S„ = 0. Then Eq. (|ic]) is used to find G n , Eq. (0) is employed 
to extract the effective medium, Eq. (||) is used to find a new local Green's function, and then Eq. (^) is used to find 
the new self-energy. The algorithm is then repeated until it converges, which usually requires only about a dozen or 
so iterations. In this contribution, we examine the half-filled case pd = J2i( n i)/N = 0.5 and Pf = J2i( w i)/N = 0.5, 
which corresponds to p = U/2 and e/ = —U/2. 

In the following section we illustrate how the Baym-Kadanoff formalism can be used to determine the dynamical 
charge susceptibility. Numerical results at half filling are presented in Section III and conclusions in Section IV. 

II. BAYM-KADANOFF FORMALISM 

The momentum-dependent susceptibility satisfies the following Dyson's equation: 

X dd (q, iuJ m , iu n \iv x ) = Xo d (^ iuj m; m)$mn 

- T Xo d (q> ium ; iui)T(iu) m , iuj n > ; ivi)x M (q, % , i^n \m), (H) 
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and the full susceptibility is found by summing over the Fermionic Matsubara frequencies X (li i v l) 
T^ mn x dd (q, iw TO ,iw n ; W;). In Eq. ( |ll|) the bare susceptibility satisfies 



n k 
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with X(q) = lim^^oo Y^t=i cosqi/d describing all of the momentum dependence of the susceptibilityBS and 

Foo (z) = [°° de^, (13) 

the Hubert transform of the nonintcracting density of states. 

The dynamical susceptibility simplifies in three cases: the chessboard case [q = (jr, it, 7r), X = —1] where 

dd,l -i \ G m -\~ G 7n -^[ . 

Xo (-l,iu m ;m) = —. — — - - ; (14) 

lLU rn + lUJ m+ l + — 2j to — L, m+ l 

the local case (X = 0) where 

Xo rf (0, iuj m ; m) = -G m G m+ i; (15) 
and the uniform case [q = (0, 0, 0), X = 1] where 

Xo {l,iui m )in) = -— — . (16) 

IVl + 2j m — l^rn+l 

We will only be interested in these three simpler cases here. Note that in instances where the denominators of Eqs. (|l4| ) 
and ( p^ ) vanish, the susceptibility is evaluated by l'Hopital's rule. 

The [-calculation of the full susceptibility requires the local irreducible vertex function. The Baym-Kadanoff ap- 
proachQu solves for the irreducible vertex function in a manner which guarantees that an approximation scheme 
maintains the conservation laws of the Hamiltonian. The procedure requires that an approximate self-energy, 
£(ri, ti, r2, T%j be given by the functional derivative of a free-energy functional, $, with respect to the full Green's 
function, G(i"i, n, r 2 , t 2 ). That is: 

E(n, n, r 2 , r 2 ) = — -, (17) 

dG(r 1 ,Ti,r 2 ,T2) 

and then differentiation of the self-energy with respect to the Green's function produces the irreducible vertex function 

r(ri,ri,r 2 ,T2;r3,T 3 ,r 3 ,T3 = -— -. (18) 

T <5G(r 3 ,T3,r3,r 3 ) 

Of course, this technique can also be used for an exact solution as we do here. 

Using the dynamical-mean-field approximation for the Falicov-Kimball model in a spatially invariant system yields 
an exact expression for the self-energy E = Gq 1 — G _1 [as seen in Eq. ([?])] and we do not need to find the appropriate 
free-energy functional Our strategy is to calculate the Green's function of the effective medium, Go(ti,t 2 ), 
when there is an additional time-dependent field, x( T )? which couples to the charge density, d(r)d^(r). [i.e. we 
add Jq &Tx{T)d^ (T)d(T) to the exponent of the evolution operator in Eq. (0)]. This time-dependent field removes 
time-translational invariance from the system, so the Green's functions now depend on t\ and r 2 separately. Using 
Go(ti,t 2 ), we evaluate the self-energy as an explicit function of the full Green's function, G(ri,r 2 ), including terms 
to linear order in x( T )- We take the derivative of the self-energy with respect to the full Green's function to obtain 
the vertex function and afterwards set the field, x( T ) to zero. The calculation differs from standard approaches, in 
that x( T ) provides a time-dependence so that Go(ti,T2) depends separately on t% and r 2 , not just on the difference, 
(ti-t 2 ). n 

To begin, we introduce an auxiliary Green's functiom (as Brandt and Urbanek did in their calculation of the 
/-spectral function) in the presence of just the charge-coupled field, Xi defined by: 
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(19) 



1 + e^eic :d- x (-)dt(-) d( -) 

where ifo = —jiSd and we expand the time-dependent field, x, in a Fourier series of the Bosonic frequencies 

X(r) = ]>>e W!T . (20) 



Noting that d(r) = e rff °<i(0)e THo = e TM <i(0), and taking into account the time ordering, we can easily perform the 
trace in Eq.(|l9|), which yields (to linear order in xi)'- 



gaux(Tl,T2] /J.) = C 
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which depends on both t\ and t 2 (not just T\ — t 2 ) because of the terms with xi=£o- Nevertheless, the following 
symmetries are readily proven: 



and 



for Tl < T 2 < Tl + (3 

g aux in + P, t 2 ; fi) = - g auxin, t 2 ; n), 



for r 2 < ti < t 2 + (3 

gauxin, n +/3;/x) = - g auxin, t 2 ;h). 



Hence we can find the Fourier transform of g a u:r in , t 2 ; /x) in terms of the Fermionic Matsubara frequencies 



dr 2 e^" ri e~ l ""' T2 ( 
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The result is that g aux iiu n , iu> n r, [j,) contains terms that are diagonal in n,n' and terms shifted off the diagonal by I 
that are proportional to xi'- 



lauxiiu n ,iu n >; (i) = 
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The Green's function for the effective medium of the atomic problem, Go(r, r'), is obtained as before by including 
a local time-dependent field, A(r — r'), that incorporates the effects of propagation on the lattice and only depends 
on the time-difference: 



Goir,T') = 



-Tr.fr,- ( e-^oJo^xWi-W^sWdin^fa) 



^o(m) 



(26) 



As the local field only contains diagonal frequency components, \(iuj n )5 nn i , the effective-medium Green's function is 
easily obtained from the auxiliary Green's function through: 



[ G 0] n ?n' = [9auxiv)] n y ~ A(iw„)<5„ 



(27) 



where the Green's function matrices are inverted in frequency coordinates, such that [GoL„/ represents the n,n' 
components of the inverse of the matrix Go(wj ra , iu> n i). A key aspect of our derivation of the vertex function, is that 
to linear order in a single frequency component of the charge-coupled field, xi, the Green's function for the effective 
medium, Go(ioj n ,iu> n '), like the auxiliary Green's function, g a uxiiu n , iu> n >; fx), only contains components in the two 
diagonals given by <!>„,„< and 5 n +i,n'- 
The full Green's function is defined by: 
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(28) 



where Z is given by Eq. (Q) and the time dependence of the d-operators is governed by H at . The Green's function 
contains a trace over /-electron states, which leads to two terms: 



-Tr d T T ( e 



G(ri,r 2 ) 



,-(3H J" drx(r)dt (r)d(r) )d t ) 
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-Tr d T T ( e 



T 2, 



(29) 



(30) 



where 7?o = —(icftd is the atomic Hamiltonian with no /-electrons, and Hi = (U — n)d'd is the atomic Hamiltonian 
for (i-electrons in the presence of one /-electron. 

We can now relate the full Green's function to the Green's function of the effective medium through the following 
matrix equation: 



G(iu> n ,iui n r) = w G (iu} n ,iu> n >) +w\ [(G x ) - U] , 



(31) 



where wq — 1 — w\ and w\ are given by Eq. ^ and equal the fraction of sites with an occupancy of zero and one 
/-electron, respectively. In the above, and all following equations, U is the only matrix which is necessarily diagonal 
in frequency space. 

Eq.([3l|) can be rearranged by multiplying on the left or right by matrices like G -1 , Gq 1 and Gq 1 — U to give the 
two following matrix equations (with indices suppressed): 



Gq 2 -(U+ G-^Go 1 + (1 - w^UG- 1 = 0, 
G o 2 - Gq\U + G- 1 ) + (1 - w^UG- 1 = 0. 



(32) 
(33) 



Adding these two equations together and collecting terms (noting the noncommutativity and hence the ordering of 
the matrices) then yields the following quadratic matrix equation: 



G^-^U + G- 1 ) 



1 2 



-\{U+ G- 1 ) 2 + (1 - w^UG- 1 = 0. 



Substitution for Go by E, using Dyson's equation in the form 

T,(iuj n ,iuj m ) = [Go l ] nm - [G _1 ] 

yields: 



(34) 



(35) 



S+-(G- 1 -[/) 



[U 2 + 2(2tei - ljUG- 1 + G- 2 ] , 



(36) 



where each term in G and £ is a matrix in frequency-space and U multiplies the identity matrix. Note that in the 
limit xi = 0j Eq.(p6[) becomes diagonal, and reduces to the conventional quadratic expression for the self-energy in 
terms of G(itu n ) and w\ as first determined by Brandt and Mielscba. 
The irreducible vertex function is defined in frequency space by: 



f-P rP+Ti pf3 

T(iuj n ,ioj m ;ioj n >,iuj m i) = j dn J dr 2 J dry 



&T 2 <e 711 " 1 m 2 r{Tl,T 2 ]Tl>,T2>). (37) 



As the vertex function is independent of the absolute time, we can change variables: T\ i— » t\ — T2 1 ,T2 > t 2 — T2> , Ty t— > 
tv — t 2 ' so that the T2 1 integral yields a delta function which requires iu m ' = i^n — iu m + iuj n > . Hence, given the 
difference between two Matsubara frequencies, iuj m — iui n — ivi is a Bosonic frequency, the vertex function becomes 



1 SY.(iuj n ,iuj n + ivi) 
T 5G(iu> n >, iui n i + ivi) 



(38) 



5 



The problem is greatly simplified by calculating £„. m to first order in x, for a single frequency component, \i- I n 
this case, both S„ jm and G„. m contain diagonal terms, and off-diagonal terms linear in xi where m = n + I, so that: 



G„ 



Gn&n,m ~t~ G n 3n-\-l,m- 



(39) 



We are interested in the. ease where I ^ 0, as the zero- frequency response (corresponding to a shift of the chemical 
potential) is well-knownB'tffl. 

Next, we calculate the diagonal and off-diagonal pieces of the self-energy. Eq.(|36|) simplifies from a full quadratic 
matrix equation, to the following coupled equations: 
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2G n G n +i 



2G, 
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2G 
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G n 



4G n G n 



+i 



2(1 - 2w 1 )U 



(G„ + G n +i) 



G n G n 



+i 



(40) 
(41) 



Eq. ( p0[ ) contains only diagonal terms with no dependence on the finite- frequency field xi > hence £ n is unchanged from 
its value when xi — 0- We concentrate on rearranging Eq.(fhij) to obtain the frequency-dependent response. First, we 
multiply both sides of Eq. ( jil] ) by 



U ^ 1 
2 2G(x 



— I E„+; — 



U 



2G n 



(42) 



and use Eq.(|40|) to replace the quadratic terms on the left, giving: 



Gn 



2G n G n 



+i 
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G n 
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E n 



1 



2G n 



— E n -U ~ 



1 



2G n 



(43) 



G n G n +l 

The above form gives rise to considerable cancellation between the two sides, leading to the amazingly simple result 



£«. — Gr, 



G n — G n +i 



(44) 



So the only terms in S n m that depend on £l finite frequency component, ii^i, of the ^ field axe S n n _[_; — ^n* 

The only 

terms which survive the limit, xi ^ after differentiation, to give the vertex function, are: 



1 5T,(iujn,iujn + i) 1 SY, n !£„-£„ 
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The .above simple result found from this Baym-Kadanoff analysis is identical to that found by diagrammatic calcu- 
lationsu which give the more complicated form: 



Y{ioj n ,iLO n ;ivi) = 



WqWiU 2 



T (1 + G„E„) [1 + G n (En U)] (1 + Gn+l^n+l) [1 + G n+l (Z n +l ~ U)] + W oWl U 2 G n G n+l 



(46) 



The proof of this equivalence is presented in the Appendix. 

Since the vertex function is diagonal, we have a simple form for the dynamical charge susceptibility 



1 (X;iv l ^0)=Tj2 



Xo(X,iu) m ;ivi) 



1 + xo{X,iu} m ;ii/i) 



G m -G„ 



In the limit q = (X = 1), substitution of Xq from Eq. (|16|) into Eq. fl47|) gives 



eld 



G m — G 

in 



m-\-l 



= o, 



(47) 



(48) 



G 



as expected by our symmetry arguments. Notice that Xo d does not vanish; it is the vertex corrections that force the 
full susceptibility to vanish. 

The final step in our formalism development is to perform the analytic continuation of Eq. (|47]) from the imaginary to 
the real axis. Our method is not completely rigorous, because we are unable to verify necessary analyticity arguments 
as described below, but we compare the direct calculation of the susceptibility on the imaginary axis to the inferred 
value there that is found by using the spectral formula from the real axis. In nearly all cases, we have accuracy to 
better than one part in 1000, which supports, a posteriori, that our technique is valid. 

We will consider only the case with X = — 1 here, since the X = 1 susceptibility is trivial and since the X = case is 
simpler and can easily be worked out by following the same steps we use for the chessboard case. We start by replacing 
the summation over Matsubara frequencies in Eq. ( [47|) by three contour integrals, which are chosen as in Figure 1 
to encircle all of the simple poles at the Fermionic Matsubara frequencies, and no other poles in the corresponding 
integrals [in other words, the only poles contributing are the simple poles from f(u>) below]. The susceptibility then 
becomes (for simplicity, we consider vi > here) 

GR{u)+GnXui+iui) 



2tt 
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2^ 



c 2 



'c 3 



1 - 
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'G A (u)+iv l ) 



(49) 



where f(uo) = 1/[1 + exp(/3w)] is the Fermi function and the subscript R or A refers to the retarded or advanced 
Green's function (or self-energy). The choices of the subscripts are so that the Green's functions and self-energies 
are analytic within regions 1, 2, or 3. Hence the contours can be deformed until they run parallel to the real axis, as 
shown in Figure 2. We are making an assumption that there are no other poles present within these regions, when we 
deform the contours. In fact, if we make an analytic continuation of just the bare susceptibility, then after continuing 
ivi to v + iS we do find a pole that lies in region 2 just below the real axis (at uj = — v/2 — iS), but the residue of this 
pole vanishes when v lies on the real axis. It is more difficult to make such an analysis for the full susceptibility, so 
we rely instead on the comparison with the direct calculation on the imaginary axis. 

When we evaluate the integrals along the lines indicated in Figure 2, we will evaluate the Fermi function at w- ivi, 
which we set equal to f{uj) before continuing the v frequency Then we can evaluate the final result, which becomes 

C G(u)+G(u>+v) 



00 I 1 2w+2/. t +iy-S(w)-S(w+iy) G(uj)~G(uj+v) 

G' ju) + G' {uj+v) 
,, , 2u+2p+i/-E«(tj)-S*(u+i') 

~ J\ W ~ V ) G'(uj)+G*{lo+v) 

1 2w+2ft+K-£*(u;)-E (w+i/) G*{w)-G*\iO+v) 

G* (uj)+G(u+v) \ 

r r/ \ r/ M 2w + 2^l + I^-S* (w)-S(w + l/) I /- a n 

-\f(u)-f(u>-v)] - G -M+G(a,+„) s-m-so^) f' ( 5 °) 

2uj+2ti+v-Y,*{u))-Y,(uj+v) G*(J)-G(u)+v) ) 

and we replaced the advanced functions by the complex conjugate of the retarded Green's function (which is valid on 
the real axis). 



III. NUMERICAL RESULTS 



We performed a number of different numerical calculations of the charge response. To begin, we summarize the 
d-electron spectral functions in the weak-coupling (U = 0.25) limit, the moderate coupling limit (U = 1) and the 
strong-coupling limit (U = 4), which are illustrated in Figure 3. Nate that the interacting density of states is 
temperature independent for this model (as first shown by Van Dongenlij). In the weak-coupling case (U = 0.25) the 
density of states looks Gaussian, with only small modifications due to the interactions. At moderate coupling, (U = 1), 
we find that a pseudogap appears in the density of states at the chemical potential. Note that this is a correlation- 
induced pseudogap and is not resulting from the charge-density-wave order of the ground state, since these results are 
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for the high-temperature homogeneous phase. Finally, a true gap develops in strong coupling (U = 4), which rapidly 
approaches U as U increases in magnitude. In Figure 4 we show plots of the static chessboard susceptibility for the 
three different values of U and inset a plot of the transition temperature as a function of U/(U + t*). The chain-dotted 
lines are guides to the eye for T c at the three values of U. The transition temperature has a classic form — increasing 
for weak-C/, reaching a maximum at moderate U and then decreasing in the strong-coupling regime. Notice how the 
static susceptibility diverges as we approach T c because we are in the thermodynamic limit. In Figure 5 we show 
the dynamical chessboard susceptibility at five different temperatures ranging from well above T c to just above T c . 
Notice how there is a temperature dependence to the charge susceptibility, but how there is no indication of the phase 
transition seen in the static susceptibility (recall the susceptibility is discontinuous at v = and we are only plotting 
the continuous piece of the susceptibility). In Figure 6 we show the same plot for the local dynamical susceptibility. 
Its behavior is quite similar to that of the chessboard susceptibility, but with a less marked temperature dependence. 
The expected peak in the imaginary part at low frequencies due to quasiparticle excitations can be seen, and it has a 
classic "triangular" shape as expected in a nonintcracting system (recall the local susceptibility is the average over all 
momentum vectors, so we expect a linear contribution at small frequency and a curved drop off at high frequencies 
as we reach the approximate half bandwidth). The case of moderate coupling is shown in Figures 7 and 8. Here, the 
chessboard susceptibility actually decreases as T is lowered because of the pseudogap in the density of states. Once 
again the temperature dependence of the local susceptibility is much less than the chessboard susceptibility (and we 
are including data below the T c to the charge-density- wave instability). The development from a single low-energy 
peak in the real part, to two nearly overlapping peaks (corresponding to charge transfer excitations), is clear even at 
moderate coupling where the density of states has only a pseudogap. Note once again that these peak developments 
arise from correlation effects and not from an underlying charge-density-wave order because these calculations are all 
performed within the homogeneous phase (even when we go below T c ). 

In Figures 9 and 10 we show the results at strong coupling U = 4. The charge susceptibility shows a number of 
interesting new features in the strong-coupling regime. The charge-transfer peaks, which are easiest seen as broad 
peaks in the imaginary part of the susceptibility, are present at the expected locations of ±U. At high temperatures, 
the system has a large peak in the real part of the susceptibility near zero frequency, which rapidly decreases as T is 
lowered, and becomes unnoticable at low temperatures. The imaginary part goes from having a linear behavior at low 
frequency to being essentially zero. Hence, in addition to the energy scale of the order of U (corresponding to the large 
charge-transfer peaks in the imaginary part of the susceptibility) there is a low-energy scale on the order of t* 2 /U 
that determines the low-temperature evolution of the charge excitations and the energy scale of those low-energy 
excitations. The high-energy features are "frozen in" at a high temperature scale, and there is strong temperature 
evolution in the low-frequency regime until T becomes smaller than the low-energy scale where the low-energy features 
are "frozen in." We believe these low-energy charge excitations are associated with virtual hopping processes, where 
an electron virtually hops onto a site occupied by a localized electron, and then hops off that site. Such processes 
turn off at temperatures below t* 2 /U. Once again the dependence of the chessboard and the local susceptibilities are 
similar and there is no signal of the phase transition in the dynamical piece of the charge susceptibility. We find that 
we lose some accuracy in our calculations at strong coupling, perhaps arising from the gap in the density of states and 
its effect on the low-temperature charge dynamics. Another possibility is that there are additional poles and resides 
that need to be taken into account in our "approximate" analytic continuation. One illustration of the numerical 
difficulties is that we no longer have an accuracy of one part in a thousand (it is typically a few parts in a thousand) 
when we evaluate the charge susceptibility on the imaginary axis by using the spectral formula and comparing with 
a direct calculation on the imaginary axis. Another illustration is that at the lowest temperatures the imaginary 
part of the spectral function becomes slightly negative for a frequency range in the vicinity of 0.5 < v < 2 (where 
the imaginary part is approximately —0.0003), which must be a numerical artifact due to the fact that the spectral 
function is known to be nonnegative for positive values of the frequency. 

IV. CONCLUSIONS 

We have presented a nontrivial example of the Baym-Kadanoff formalism to derive the dynamical susceptibilities 
of the Falicov-Kimball model. The dynamical susceptibility breaks into two pieces, a static piece that reflects the 
coupling of the system to the charge-density-wave distortions, and a frequency-dependent piece, which is less affected 
by any underlying charge-density-wave order. We developed the formalism to evaluate this frequency-dependent piece 
of the susceptibility exactly on the imaginary axis. We also developed an analytic-continuation scheme that we were 
unable to establish with full rigor due to the possibility of a small contribution from neglected residues of poles that 
could arise within the analytic-continuation procedure, but which is nevertheless quite accurate over a wide range of 
parameters. 
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We find a number of interesting properties of this model: (1) the static and dynamical susceptibilities are 
decoupled — the static susceptibility diverges at T c , but no signal of this phase transition is seen in the dynamical 
susceptibility; (2) the momentum dependence of the dynamical susceptibility is not too strong around the chessboard 
point — we see little variation between the X = — 1 and X — susceptibilities, but the X — 1 susceptibility vanishes by 
symmetry, so momentum dependence is stronger near the Brillouin zone center; and (3) the dynamical susceptibilities 
do possess temperature dependence even though the interacting density of states is temperature independent. The 
temperature dependence is most striking in the strong-coupling regime where a low-energy scale on the order of t* 2 /U 
arises and the charge susceptibility has a strong temperature dependence near this approximate temperature. 

Our application of the Baym-Kadanoff formalism provides an interesting counterpoint to the diagrammatic deriva- 
tion about the atomic limit, and can be applied to other related problems where the so-called "static approximation" 
is exact. We have not examined any effects away from half-filling or at incommensurate wavevectors here. We leave 
those tasks to another publication. 
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APPENDIX A: EQUIVALENCE TO THE ATOMIC DIAGRAMMATIC APPROACH 

To compare our result with the diagrammatic approach^, we must begin with Shvaika's result for the frequency- 
dependent piece of the reducible vertex function, 



Y(iu) n ,iuj n i\ivi ^ 0) = 

1 <W U 2 wi{l-wi) 



T G n G n >G n+ iG n > + i (iu) n + fj, - A n ) (iw n + fi-X n -U) (iuj n+ i +fi- X n +i) {i^n+i + M ~ K+i ~ U) ' 
By making the replacement, iu) n + /x — A„ = [Go(iu; n )] 1 = (G n ) 1 + S n , Eq. flAip is equivalent to: 
1 S nn >w wiU 2 



(Al) 



(A2) 



T(l + G n £„)[l + G n (5jt» 

U)} (1 + G„ +i S„ +i ) [1 + G n+l {H n+l - U)] " 
Solving Dyson's equation for the irreducible vertex function, T(iu n , iu) n i\ ivi ^ 0): 

T{iu} n ,iu) n i;ivi ^ 0) = T(iw n ,iw n i;ivi ^ 0) - T ^T) T(iw n , ico n " ; wi ^ 0)Gn"G n i> + ir(ioj n »,i(j n r,iv l ^ 0), (A3) 

n" 

then yields Eq. @. 

In order to prove that our simple result, Eq.(|4^) and the diagrammatic result, Eq.(^6|) are equivalent to each other, 
we first demonstrate the existence of a useful identity. Eq.(^o|) can be expanded and multiplied by G n to yield: 

Wl U = E n (l-UG n + E n G n ), (A4) 
= S n +i (l — U G n +i + Y> n +lGn+l) , (A5) 



where the second equation is evaluated equivalently, with bh n + I. Eq. (A4) is used twice to simplify the product 
given below: 



(1 + G„S n ) [1 + G„(E„ - U)] = 1 + G„E„ - UG n + w x U, (A6) 



= wxU[^- + G n ), (A7) 



and similarly for n i— * n + I . 

Hence the expression on the right of Eq. (^) is rewritten as: 

1 WqWiU 2 



w l U2 (sr + G ») (sin + Gn+l ) + GnG n+ iw oWl U 2 



(A8) 



and multiplication of top and bottom by (E„ — £„+;)£„£„+; leads to: 
1 (E n - E„+;)E n E n+ ;(l - wi) 



T wi [E n — E n +i + G n E n (E„ — S n +i) + G ri +iE n +;(E ri — E„+i)] + G„G n +zE ra E ra +i(E n — E„+;) 



(A9) 



Now Eqs. JA4HA5| ) are further used to make the replacements G rl E^ i— * wiU — E n + C/G rl E„ and G n +iE^_|_ ; 
w\U — E n +; + UG n +i^n+i for all quadratic terms in the denominator. Numerous cancellations then lead to: 



T{iu) n ,iw n ;ivi ^ 0) = 



1 



(E„ — E n+ i)E„E n+ ;(l — wi) 



T wi (E n +iE n G r j+i — E n -|_;E ra G n ) - En+/E n Gn+^ ~t~ E n +/E n Gn 

1 E„ — E n _|_; 



T G n ~~ G n 



+1 



(A10) 
(All) 



which yields our final result, Eq. (p 
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Freericks and Miller, Phys Rev B, Figure 1 

FIG. 1. Contour integral for evaluating the Matsubara frequency summations of the charge susceptibility. The x's mark 
the locations of Fermionic Matsubara frequencies. The contours enclose all Matsubara frequencies, but no other poles of the 
system. Note that we divide the complex plane into three regions: (i) region 1, where the imaginary part is greater than zero; 
(ii) region 2, where the imaginary part lies between zero and —ivi\ and (iii) region 3, where the imaginary part is less than 
— ivi- 



10 



X 



ys Region 1 

X 



-,- CO 

X 

X 

Region 2 

X 

X ^ 

" CO-iV| 

X 
^/ 

Region 3 

X 
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FIG. 2. Deformation of the contours needed for evaluation of the susceptibility on the real axis. The integrations are parallel 
to the real axis. 




Frequency (cj/t*) 

Freericks and Miller, Phys. Rev. B, Fig. 3 

FIG. 3. Interacting density of states for the spinless Falicov-Kimball model at half filling and U = 0.25 (solid line), U = 1 
(dotted line), (7 = 4 (dashed line). Note how the density of states evolves from being essentially Gaussian to developing a 
pseudogap and then a real gap as the interaction strength is increased. The density of states is temperature independent for 
this model. 
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Freericks and Miller, Phys. Rev. B, Fig. 4 

FIG. 4. Static charge susceptibility for the chessboard charge-density-wave plotted versus temperature for three values of U : 
(i) U = 0.25, solid line, T c = 0.0153; (ii) U = 1, dashed line, T c = 0.0608; and (iii) U = 4, dotted line, T c = 0.0547. Inset is the 
charge-density-wave transition temperature plotted as a function of the interaction strength U/(U + t*). 
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Freericks and Miller, Phys. Rev. B, Fig. 5 

FIG. 5. Dynamical chessboard charge- density- wave susceptibility as a function of temperature for the case U = 0.25. The 
real part is plotted in (a) and the imaginary part in (b). Five different temperatures are shown: T = 0.5 (solid line), T = 0.25 
(dotted line), T = 0.125 (dashed line), T = 0.0625 (chain-dotted line) and T = 0.03125 (chain-double-dotted line). The real part 
of the susceptibility becomes more sharply peaked at zero frequency, but shows no sign of diverging (because of the decoupling 
of the static and dynamic susceptibilities) . The imaginary part sharpens as T is lowered, and has a shape reminiscent of that 
of a nonintcracting system arising from the quasiparticle excitations. 



12 



o" 

II 

X 

y< 

CD 



0.5 



o = 



^ °- 5 
o" 

11 o 

X u 

E-0.5 



1 1 1 1 

A 


1 1 


' '(a)" 








/ ^ 








1 


^— — — - 


f ■' 1 ' ' ' 1 


' 1 


■■ 'i 






T=0.5t* 

T=0.25t* 

T=0.125t 

T=0.0625t 

T=0.03125t 

I 


■ X/ 

1,1,1 




i 



-4 -2 2 4 

Frequency (v/t*) 

Freericks and Miller, Phys. Rev. B, Fiq. 6 

FIG. 6. Dynamical local charge- density- wave susceptibility as a function of temperature for the case U = 0.25. The real part 
is plotted in (a) and the imaginary part in (b). Five different temperatures are shown: T = 0.5 (solid line), T = 0.25 (dotted 
line), T = 0.125 (dashed line), T = 0.0625 (chain-dotted line) and T = 0.03125 (chain-double-dotted line). The form of the 
susceptibility is similar to the chessboard case, except the system reaches the low-temperature limit much more rapidly, and 
the imaginary part is pushed somewhat farther from zero frequency. 
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FIG. 7. Dynamical chessboard charge-density- wave susceptibility as a function of temperature for the case U = 1. The real 
part is plotted in (a) and the imaginary part in (b). Five different temperatures are shown: T = 0.5 (solid line), T = 0.25 
(dotted line), T = 0.125 (dashed line), T = 0.0625 (chain-dotted line) and T = 0.03125 (chain-double-dotted line). The real 
part of the susceptibility initially increases as T is lowered at low frequencies, but then decreases and starts to develop additional 
peaks as the effects of the pseudogap are felt. 
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FIG. 8. Dynamical local charge-density-wave susceptibility as a function of temperature for the case U = 1. The real part 
is plotted in (a) and the imaginary part in (b). Five different temperatures are shown: T = 0.5 (solid line), T = 0.25 (dotted 
line), T — 0.125 (dashed line), T — 0.0625 (chain-dotted line) and T = 0.03125 (chain-double-dotted line). We once again see 
a similarity with the chessboard phase, but the temperature dependence is much reduced here. 
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Freericks and Miller, Phys. Rev. B, Fiq. 9 
FIG. 9. Dynamical chessboard charge-density-wave susceptibility as a function of temperature for the case U = 4. The real 
part is plotted in (a) and the imaginary part in (b). Five different temperatures are shown: T = 0.5 (solid line), T = 0.25 
(dotted line), T = 0.125 (dashed line), T = 0.0625 (chain-dotted line) and T = 0.03125 (chain-double-dotted line). Note how a 
substantial low-energy feature, present at T = 0.5 disappears rapidly as T is lowered. The temperature scale for this evolution 
is on the order of T ^ t* 2 /U = 0.25. For temperatures below this energy scale the system's behavior is governed most strongly 
by the energy gap in the density of states. We find that the imaginary part of the charge susceptibility goes slightly negative 
over a frequency range around 0.5 < v < 2, which we believe to be an artifact of the numerical analytic continuation procedure 
in the presence of a gap. 
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FIG. 10. Dynamical local charge-density-wave susceptibility as a function of temperature for the case U = 4. The real part 
is plotted in (a) and the imaginary part in (b). Five different temperatures are shown: T — 0.5 (solid line), T = 0.25 (dotted 
line), T = 0.125 (dashed line), T = 0.0625 (chain-dotted line) and T = 0.03125 (chain-double-dotted line). The results here 
are similar to those for the chessboard susceptibility, except here we have a slower temperature dependence, asav = peak 
in the real part of the charge susceptibility can still be seen at T — 0.25, and it has already vanished in the chessboard case. 
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